Boundary and expansion effects on two-pion correlation functions in relativistic 

heavy-ion collisions 
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We examine the effects that a confining boundary together with hydrodynamical expansion play 
on two-pion distributions in relativistic heavy-ion collisions. We show that the effects arise from the 
introduction of further correlations due both to collective motion and the system's finite size. As is 
well known, the former leads to a reduction in the apparent source radius with increasing average 
pair momentum K. However, for small K, the presence of the boundary leads to a decrease of the 
apparent source radius with decreasing K. These two competing effects produce a maximum for 
the effective source radius as a function of K. 

PACS numbers: 25.75.-q, 25.75.Gz, 25.75.Ld 



I. INTRODUCTION 

In recent years, much experimental effort has been de- 
voted to the production of a state of matter under the ex- 
treme conditions of high baryonic densities and/or tem- 
peratures in relativistic heavy-ion collisions. The main 
drive for such effort is the expectation to produce the so 
called quark-gluon plasma (QGP) where the fundamental 
QCD degrees of freedom are not confined within a single 
nucleon but rather over a larger volume of the order of 
the dimensions of a nucleus. While the properties of the 
QGP have been the subject of intense theoretical study 
and debate, much less attention has been paid to the 
hadronization process, or to the properties of hadronic 
matter at high temperatures and densities. An under- 
standing of these is needed for a correct interpretation 
of the signals originating from the different stages of the 
collision, both for a clear distinction of a possible QGP 
formation but also as an interesting subject of study on 
its own. 

The most abundantly produced hadrons in relativistic 
heavy- ion collisions arc pions. Typically, the number of 
pions produced one unit around central rapidity in cen- 
tral Au+Au collisions at energies of order 10A GeV is 
dN^/dy ~ 300 pj. Under the assumption that the trans- 
verse dimensions of the system formed at freeze out are 
of order of the transverse size of an Au nucleus and that 
the typical pion formation time is of order 1 fm, this 
multiplicity implies that the average pion separation d at 
freeze out in the central rapidity region is of order d ~ 0.6 
fm which is less than the average range of the pion strong 
interaction, d s ~ 1.4 fm. 

Some of the possible consequences of this large pion 
density produced in relativistic heavy-ion collisions where 
first studied by Shuryak |2j who coined the term pion liq- 
uid to refer to the situation where the pion system could 
not be thought of as existing as a hadron gas but rather, 
that its properties resembled more those of a liquid of 
quasipions. In particular, as one of the main character- 



istics of liquids is the appearance of a surface tension, 
such state of hadron matter could give rise to a confining 
boundary that acted as a reflecting surface that could af- 
fect thepion distributions just before freeze out. More re- 
cently H , it has been proposed that the equation of state 
of pion matter could give rise to a phase transition from 
a gas phase to a more dense phase as the temperature 
rises close to the temperatures expected to be achieved 
in relativistic heavy-ion collisions, thus introducing the 
concept of a hot pion liquid. 

From a phenomenological point of view and disregard- 
ing the details of the reflection processes (which presum- 
ably depend on the energy of the incident particle), the 
development of a surface tension can be incorporated by 
imposing a sharp boundary for the pion wave functions 
to evolve just before freeze out. As a consequence of the 
finite size of the system during this stage of its evolution, 
the energy states form a discrete set. 

An important difference between statistical systems 
with and without a boundary is the different density of 
states at increasing energies, being larger in the case of 
the former, as illustrated in Fig. 1. The density of states 
of a finite system approaches that of an unbound one as 
the size of the system is increased. The above character- 
istic implies that the transverse inclusive spectrum cal- 
culated within a boundary model will exhibit a concave 
shape at high transverse momentum and could poten- 
tially explain the increase of the pion transverse distri- 
bution |i) in this kinematical region The increase of 
the pion distribution at low transverse momentum could 
also be explained within the same context by considering 
the finite chemical potential associated with the mean 
pion multiplicity per event p] . 

More recently, it was realized that an important miss- 
ing ingredient in the description of the transverse pion 
spectra within a boundary model was the proper inclu- 
sion of hydrodynamical expansion. The phenomenologi- 
cal description in terms of a bound, expanding pion sys- 
tem has been named the expanding pion liquid model 
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and was developed in Refs. |7) and successfully applied to 
the description of the experimental mid-rapidity, trans- 
verse pion spectra in central Au+Au collisions at 11. QA 
GeV/c [Q. A further natural test ground for the model 
is the study of multiparticle correlations, in particular 
two-pion correlations. 
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FIG. 1. Number of states for a spherically symmetric sys- 
tem with a sharp boundary at R — 6 fm (solid line) compared 
to a system without boundary (dotted line). Notice how the 
presence of the boundary makes the number of states grow 
faster at large energies. 

A step in this direction has been taken in Ref. || where 
the effects of a pure boundary model (i.e. without collec- 
tive expansion) have been introduced in the description 
of the two-pion correlation function. In this work, we 
incorporate also the effects of hydrodynamical expansion 
in the calculation of the two-pion correlation function. 

It is well known that in the study of two-particle cor- 
relation functions, the effective size of a system without 
a boundary decreases as the average pair momentum is 
increased || when considering the effects of hydrodynam- 
ical expansion. Physically, this effect is due to the fact 
that as the average pair momentum increases, the parti- 
cles in the pair are more likely to be emitted from points 
close in space. This can also be regarded as the introduc- 
tion of a further correlation in phase space for the emitted 
particles that effectively reduces the size of the emission 
region |l(J] ■ On the other hand, as emphasized in Ref. || , 
for emission volumes of order of a few average pair wave- 
lengths it is imperative to consider a full quantum density 
matrix in the description of particle distributions since as 
the average pair momentum decreases, the correlation re- 
gion as a function of relative momentum increases thus 
effectively reducing the apparent size of the system. As 
we will show in this work, these two competing effects 
produce a maximum in the effective size of the system at 
a finite value of the average pair momentum. 

This paper is organized as follows: Section || is devoted 
to a recollection of the description of a bound, expand- 
ing pion system at freeze out. In Sec. |ffl|, we write the 
two-pion correlation function in terms of the discrete set 



perform a systematic analysis of the two-pion correlation 
function in terms of the different parameters involved. 
We pay particular attention to the behavior of the ef- 
fective system's radius as a function of the average pair 
momentum comparing the results to what would be ex- 
pected in the case of an expanding and unbound system 
and a pure boundary model. Finally we conclude and 
discuss our results in Sec. [v|. 



II. THE EXPANDING PION LIQUID 

When the system of pions can be considered as con- 
fined, its wave functions satisfy a given condition a the 
boundary. In order to compare the results with the ob- 
served particle distributions, the shape of the assumed 
confining volume could play an important role. It has 
long been known that the particle momentum distribu- 
tions are somewhat forward-backward peaked, particu- 
larly at energies of the SPS pH , even in the case of cen- 
tral collisions. Nevertheless, for the sake of simplicity 
and concreteness, here we will assume that the confining 
volume is spherical. Comparison of the model with data 
will become better in the central rapidity region where an 
asymmetry between longitudinal and transverse expan- 
sion is less important than in the fragmentation region. 
Let us emphasize that an spherically symmetric model is 
not essential to the basic physics discussed here and can 
thus be relaxed at the expense of additional computing 
time. This could however be necessary when comparing 
to data away from mid-rapidity. 

To incorporate the effects of an hydrodynamical flow, 
we notice that this ordered motion can be represented 
by a four-velocity field u M = 7(r)[l, v(r)], where 7 is the 
Lorentz-gamma factor and v(r) is the velocity vector. 
This field represents a redistribution of momentum in 
each of the fluid cells, as viewed from a given reference 
frame (the center of mass in this case), becoming centered 
around the momentum associated with the velocity of 
the fluid element. This behavior can be described by the 
substitution p^ — ► p^ — mu^ , where m is the pion mass. 
The term mu^ represents the collective momentum of the 
given pion fluid element. 

The eigenf unctions of the confined, expanding system 
of pions are thus obtained as the solutions to the equation 
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(1) 



where we look for the stationary states subject to the 
boundary condition 



ff>(\T\=R,t)=0. 



(2) 



of eigenfunctions for this kind of systems. In Sec. IV we 



We consider a parametrization of the velocity vector v(r) 
that scales with the distance from the center of the fire- 
ball 
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v(r) = , 



(3) 



where the parameter < /3 < 1 represents the surface 
fireball hydrodynamical velocity. Correspondingly, the 
explicit expression for 7(7") becomes 



7(r) 



1 



y/l - (3 2 r 2 /R 2 



(4) 



Since Eq. ([y), with the gamma factor given by Eq. (Q), 
can only be solved numerically, we resort to approximate 
the function 7 by the first terms of its Taylor expansion 



7(r) ~ 1 + 



2 i? 2 ' 



(5) 



which is valid for not too large values of (3. With this 
approximation, Eq. ([!]) becomes an equation for particles 
moving in a spherical harmonic well with rigid boundaries 
and can be solved analytically. The stationary states are 

^ nlm , (j,, t) = _^ =e -<^Wr'/(ja)ym' (f ^1^/2^ 

xifi (£±y*> - + 3/2; a 2 ,r 2 ), (6) 

where liq is a confluent hypergeometric function and 
is a spherical harmonic. The quantities A n i are the 
normalization constants and are found from the condition 
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(7) 



The parameters a n i and e n i are related to the energy 
eigenvalues E n i by 



o£, = m(E nl - m)/3 2 /R 2 
4; = E nl (E nl - 2m) , 

with i?„; given as the solutions to 



(i + 3/2) 



4a 2 



■,Z + 3/2;a^ =0 



(8) 



(9) 



Equation (||) , along with the energy eigenvalues and def- 
initions in Eqs. (@), (|J), constitute the set of (properly 
normalized) cigenfunctions in terms of which the vari- 
ous multi-particle distributions can be expressed. The 
system's finite size and the strength of the collective ex- 
pansion are given in terms of the parameters R and /3, 
respectively. In order to further proceed, it is necessary 
to specify the kind of ensemble that describes the statis- 
tical properties of the pion system. 



III. TWO-PION CORRELATION FUNCTION 

In order to describe the situation where equilibrium 
has been attained (which we assume here), the proper 



statistical distribution for the ensemble is thermal. For 
the purposes of this section, we closely follow Ref. || to 
where we refer the reader for details. Let A represent 
the set of quantum numbers {nlm'}. The corresponding 
occupation number N\ for a given state is written as 
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exp(£ A - fx)/T - 1 



(10) 



where T is the system's temperature and fi the chemical 
potential, related to the average total multiplicity N by 



i 



exp(£ A - y)/T - 1 



(11) 



Let ip\(p) represent the Fourier transformed wave func- 
tion for the state with quantum numbers A, namely 



Va(p) 



d 3 r e- ip r ^ x (r) . 



With the normalization adopted in Eq. (|6[ 
momentum distribution can be written as 



(12) 

, the one-pion 



Pi(p) 



d 3 N 



d 3 p 

jl^wi(p)^(p) 



(13) 



Similarly, and under the assumption of a complete factor- 
ization of the two-particle density matrix, the two pion 
momentum distribution can be written as 



-P2(Pl,P2) 



d 6 N 
d 3 pid 3 p 2 
Pi(pi)Pi(p 2 ) 

7 -i-^2 J B A JV A ^(pi)VA(P2; 



(14) 



from where the two-pion correlation function Gi can be 
written, in terms of P\ and Pi , as 



C 2 (pi,p 2 ) : 



P2(Pl,P2) 



Pl(pi)Pl(p 2 ) 

=1+ * 

Y,ExN x |Va(Pi)| 2 5>a^Va \MP2)\ 2 

A A 

(15) 

Notice that as a consequence of the factorization assump- 
tion, the correlation function is such that C*2(p,p) = 2. 
This is usually referred to as the completely chaotic pion 
production scenario |12|, which is the situation expected 
to occur in a heavy ion collision, given the considerable 
rescattering experienced by pions in the production re- 
gion. In contrast, if the particles were produced com- 
pletely coherently, they would occupy a pure quantum 
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state and the two-pion momentum distribution would be 
simply the product of two single-pion momentum distri- 
butions, leading to the absence of the HBT effect. 
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FIG. 2. (a) : C 2 (q) for a fixed value of K = 500 MeV and for 
R — 6 fm (solid line) and R — 10 fm (dotted line). The width 
of C2 decreases as R is increased, (b): C2{q) for a fixed value 
of R = 6 fin and for K = 300 MeV (solid line), K = 500 MeV 
(dotted line) and K = 700 MeV (thick solid line). For the 
chosen values of K the width of C2 increases as K is increased, 
(c): 62(9) for a fixed value of R — 6 fm and for K — 50 MeV 
(solid line), K = 100 MeV (dotted line) and K = 200 MeV 
(thick solid line). For the chosen values of K the width of 
C2 decreases as K is increased. In all cases, the temperature, 
surface expansion velocity and chemical potential have been 
held fixed to T = 120 MeV, /3 = 0.5c and fx = respectively. 



IV. THE EFFECTIVE RADIUS 

Armed with the eigenfunctions describing the confined 
and expanding pion system, Eq. (0) and with the ex- 
plicit expression for the two-pion correlation function in 
Eq. ( [l5| ) , it is possible to perform an analysis to describe 
the behavior of C2 as a function of the several variables 
and parameters involved. For the spherically symmetric 
problem described here, the correlation function depends 
on the angle between the pion pair momenta. For the 
sake of simplicity, let us consider the case in which both 
momenta pi and P2 are parallel. In this case, the sum- 
mation in the numerator of the second term in Eq. jl5| ) 
can be simplified, with the aid of the addition theorem 
for the spherical harmonics 

1 2/4-1 
Y, Yr(pi)Yr*(fr) = -±-. (16) 

m'= — I 

Figure 2 shows the behavior of Ci(q) as a function 
of q, the magnitude of the pair momentum difference 
q = p 2 Pi- Figure 2a shows C 2 (q) for a fixed value 
of the magnitude K = 500 MeV of the average pair mo- 
mentum K = (K2 + Ki)/2 for two values of the system 
radius R — 6, 10 fm. Notice that for a fixed K the width 
of the correlation function decreases as R is increased. 
Figure 2b shows C2 (q) for a fixed value of R = 6 fm and 
three values of K = 300, 500, 700 MeV. For the cho- 
sen values of K, the width of the correlation function 
increases as K is increased. Figure 2c shows C2 (q) for a 
fixed value of R = 6 fm and three values of K = 50, 100, 
200 MeV. Notice that in this case, the width of the cor- 
relation function decreases as K is increased. In all these 
figures, the temperature and surface expansion velocity 
have been held fixed to T = 120 MeV and (3 = 0.5c. 
These values for T and f3 are chosen in accordance to 
the analysis in Ref. [Jl3| where a correlation between the 
transverse flow velocity and the freeze out temperature 
is found in such a way that higher temperatures imply 
lower expansion velocities and vice versa. Since, at least 
for AGS energies, not too high temperatures are reached 
during the collision, the above value for T implies that for 
f3, thus, the free parameters for the model can be taken 
either as R and T or R and (3. 

It is also worth mentioning that in both of the above 
figures, the value of the chemical potential /1 appearing in 
Eq. ( |l5| ) has been fixed to pb = 0. The behavior of C2(q) 
for different values of the chemical potential is shown in 
Fig. 3. Notice that varying the chemical potential up 
to values below the onset of Bose-Einstein condensation 
(BEC) |6) does not introduce changes in the shape of the 
correlation function. This can be understood by noticing 
that even when we increase the system's density and thus 
the value of /i, we are not introducing any further cor- 
relation among the bound but otherwise non-interacting 
set of particles. This situation changes when, for a given 
temperature and system's size, /x is beyond the value to 
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allow for the ground state to accommodate a significant 
fraction of the particle population § . This is also shown 
in Fig. 3. In this case, the ground state population has 
to be treated separately from the pions coming from the 
excited states, since the former originates from a pure 
quantum state, as opposed to the assumption leading 
to Eq. (|l4|). The weight assigned to the pions coming 
from the ground state is equal to the ratio of the ground 
state population N g to the total number of pions of the 
given species N. Correspondingly, the correlation func- 
tion 6*2(9) becomes flatter and the intercept with the 
vertical axis occurs for values smaller than 2, that is, 
A = 1 — Ng/N. Thus, as the density increases, this be- 
havior signals that for the given temperature and volume, 
pions are predominantly emitted from the ground state. 
This is in agreement with the analysis in Ref. E3 . 
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FIG. 3. C2 {q) for a fixed value R = 6 fm and different val- 
ues of the ground state population fraction. For fi = (solid 
line) the ground state population is negligible and A ~ 1. 
However, when the ground state population becomes a signif- 
icant fraction of the total multiplicity, e.g. for /1 = 315 MeV, 
iV 9 /iV = 0.37 and A ~ 0.62 (dotted line). When almost all of 
the particles occupy the ground state, e.g. for n — 317.5 MeV, 
Ng/N = 0.77, A ~ 0.23 (thick solid line) and the correlation 
function becomes flatter. 



Another property of the bound and expanding system 
of pions that can be extracted from the correlation func- 
tion 6*2(3) i s the behavior of the system's effective radius 
i? e g- as a function of K. The relevant quantity to pay 
attention to is the ratio 77 = T/j(R)m/3 of the energy 
scale associated with random motion, i.e. T, to the en- 
ergy scale associated with ordered motion, i.e. ^{R)m(3. 
For K small compared to r/T -that is, when the average 
pair momentum is mostly due to random motion- i? e fj 
is an increasing function of K. This can be understood 
by noticing || that increasing K corresponds to increas- 
ing both of the momenta in the pion pair. Correspond- 
ingly, the quantum states that contribute to the momen- 
tum distributions P\ and P2 are those with increasingly 
larger quantum numbers. But, according to Eq. (10), 



these states are suppressed by their statistical weight and 
therefore, only those other states with smaller quantum 



numbers -and consequently with a larger spread in co- 
ordinate space- can contribute significantly to the corre- 
lation function, which in turn drops faster as a function 
of q with increasing K , leading to an increase in size of 
the apparent region of particle emission. However, for 
K > r/T, the collective motion dominates over the ther- 
mal component in K and the relevant physical effect that 
dictates the behavior of R e g is the correlation between 
the spatial region of emission of pions and the pair mo- 
mentum introduced by the collective expansion, in such a 
way that faster pions are more likely to be emitted from 
points close in space || leading to a reduction in size of 
the apparent region of particle emission. 
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FIG. 4. -R e fj for a fixed value R — 6 fm as a function of 
K for an expanding and bound system (solid line), an ex- 
panding and unbound system (dotted line) and a bound and 
non-expanding system (thick solid line). Notice that in the 
first case, R e ff reaches a maximum at a value of K ~ r/T, 
whereas for the second case, R e g decreases with increasing K 
and in the third case ^ e fj grows with increasing K. Also, for 
the expanding and unbound system, the parameter i? gauBa has 
been chosen in such a way as to give the same r.m.s radius 
than a rigid sphere with R — 6 fm. 

The behavior of i? e g- as a function of K is shown in 
Fig. 4 (solid line). Notice that the curve shows a maxi- 
mum for a value of K ~ r)T. The curve is obtained by 
fitting the correlation functions Ci (q) to Gaussians of the 
form 



g(q) = 1 + cxp{-q 2 Rl ff ) 



(17) 



Equation ( |l7| ) is a good description for correlation func- 
tions with large K . For small values of K, the fit is 
not as good. For comparison, also shown in Fig. 4 is 
the behavior of i? c fj for an expanding system without a 
boundary (dotted line) and for a bound system without 
the effects hydrodynamical expansion (thick solid line). 
For the former we choose a spherically symmetric phase 
space Gaussian distribution given by 



G(x,p) = e 



-7(x)(E p -vp)/T 



(18) 



with E p = -\/p 2 + m 2 and v and j(x) given by Eqs. (||) 
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and (|5|), respectively. The correlation function is given 
in terms of G(x, p) by || 



C 2 (pi,p 2 ) = l- 



\Jd 3 xG{x, Pi+E2) e -'(pi-P2) x| 2 
(f<PxG(x, Pl )) (Jd*xG(x,p 2 )) 



(19) 



For the bound and non-expanding system, the eigenfunc- 
tions are given in terms of Bessel functions of the first 
kind H (see also Ref. ||). The corresponding expression 
for R e q is obtained from that of OAq) by also fitting 
Gaussians of the form given by Eq. ( |17| ) . Notice that the 
curve representing the effective radius as a function of K 
for a bound but non-expanding system grows with K, in 
agreement with the analysis of Ref. In contrast, the 
curve representing the effective radius for an unbound 
but expanding system decreases monotonically as K is 
increased, also in agreement with the analysis of Ref. M . 



V. CONCLUSIONS 

In this work, we have studied the effects that a confin- 
ing boundary together with hydrodynamical expansion at 
freeze out, play on the two-pion correlation function, in 
the context of relativistic heavy-ion collisions. We have 
argued that the confining boundary could be produced 
as a consequence of the high pion density that can be 
achieved at freeze out in central collisions. 

We have shown that for a given system's volume and 
temperature, varying the multiplicity, and therefore the 
chemical potential, does not introduce any changes in the 
correlation function when fi is below the values for BEC. 
However, the intercept of the function 62(9) occurs for 
values less than 2 when the chemical potential is beyond 
the value to allow for BEC. A similar behavior can be ex- 
pected for a given pion density if the freeze out tempera- 
ture is below the critical temperature for BEC. However, 
this is a less likely scenario in this kind of collisions. 

We have found the behavior of C2 (q) when varying ei- 
ther R or K keeping the other variable fixed. Since the 
importance of correlation analyses rests basically on the 
information that it can provide about the physical size 
of the system produced during the collision, a main re- 
sult of the present work is the functional dependence of 
the effective system radius R e g with the magnitude of 
the average pair momentum K . We have shown that the 
interplay of the energy scales associated with collective 
and random motion, r y(R)m[3 and T, respectively, pro- 
duce a maximum for R e g at a value K ~ 77T, where 
r) = T /^{R)m(3. The physical origins of this behavior 
are the combined effects of the confining boundary and 
hydrodynamical expansion. In the regime where K is ba- 
sically due to random motion, the boundary effects are 
the most important and R e g grows as a function of K. 
However, in the regime where K is basically due to collec- 
tive expansion, the effective size of the system is dictated 
by the correlation between the points of emission and the 
pair momentum and R e g decreases as a function of K. 
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